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LARGE-DEVIATION PRINCIPLES FOR CONNECTABLE RECEIVERS IN 

WIRELESS NETWORKS 


CHRISTIAN HIRSCH, BENEDIKT JAHNEL, PAUL KEELER, AND ROBERT PATTERSON 

Abstract. We study large-deviation principles for a model of wireless networks consisting 
of Poisson point processes of transmitters and receivers, respectively. To each transmitter 
we associate a family of connectable receivers whose signal-to-interference-and-noise ratio is 
larger than a certain connectivity threshold. First, we show a large-deviation principle for the 
empirical measure of connectable receivers associated with transmitters in large boxes. Second, 
making use of the observation that the receivers connectable to the origin form a Cox point 
process, we derive a large-deviation principle for the rescaled process of these receivers as the 
connection threshold tends to zero. Finally, we show how these results can be used to develop 
importance-sampling algorithms that substantially reduce the variance for the estimation of 
probabilities of certain rare events such as users being unable to connect. 


1. Model description and main results 

We consider a stochastic-geometry model for a wireless network consisting of a family of 
transmitters and a family of receivers. Transmitters and receivers are modeled by indepen¬ 
dent homogeneous Poisson point processes X and Y in whose intensities are assumed to 
be non-zero and finite and will be denoted by Ay and Ar, respectively. For instance, we may 
think of transmitters and receivers as users participating in a device-to-device communication 
where messages need not be routed via a base station. It is believed that this form of com¬ 
munication will be a central concept in next-generation wireless networks [7]. The most basic 
requirement in the design of such networks is to guarantee satisfactory quality of service on 
average. Additionally, it is desirable to control and quantify the probability of low quality of 
service to occur. This necessitates a more detailed probabilistic analysis and the theory of large 
deviations provides the appropriate tools. 

Let us now describe the communication model. In order to determine the connection quality of 
messages sent out from a transmitter located at x E to a receiver located at y E M d , the signal- 
to-interference-and-noise ratio (SINR) has been identified to be of fundamental importance [3j. 
More precisely, we assume that signals are transmitted with some positive powers P x and decay 
according to the path-loss function £(\x — y |), where | • | denotes the Euclidean norm in 
and £ : [0, oo) -A [0, oo) is a decreasing function satisfying £{r) E o(r~ a ) for some a > d. 
In particular £ is bounded and f £{\x\)dx < oo. In addition to the deterministic decay over 
distance, the signal strength is also influenced by random fading effects that are encoded in a 
positive random variable F x y . Such fading effects can for example come from large obstacles in 
the environment or multi-path interference due to moving reflectors @J Chapter 22]. 

Furthermore, considering a signal sent out from X j, the strength of the interference experi¬ 
enced at a location y E M. d is assumed to be of the form 

I(Xi, y) = I(Xi, y, X) = w + Y, P Xj F Xj J(\Xj - y\). 
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In words, the interference strength at a given location y consists of a contribution from the 
thermal noise w > 0 and the aggregated signal strengths coming from all other transmitters. 
For notational convenience, we differ from the common convention [3] and include the thermal 
noise w in the interference term. Hence, the SINR for the transmitter Xi E X and the possible 
receiver location y E is defined as the ratio of the signal strength by the interference, i.e., 

S\m(X h y) = S\m(Xi,y,X) = Px * Fx *’V^ X i - vti . 

I{Xi,y) 

We assume that a connection can be established between Xi € X and Yj G Y if SINR(JQ, Yj) > t 
for some fixed connectivity threshold t. The importance of the SINR stems from Shannon’s law 
in information theory, which provides an explicit formula expressing the maximum possible data 
throughput in terms of SINR, see |H Chapter 16]. 



Figure 1 . Realization of the network model. Transmitters (red) are connected 
to receivers (blue) by black lines when a fixed SINR-threshold is exceeded 


In the present paper, we analyze how connectivity properties of the SINR-based network model 
described above behave in certain asymptotic regimes. First, we associate to each transmitter 
Xi the family of receivers Y^' that are connectable to Xi, i.e., 

= {Yj € Y : S\m(Xi,Yj) > t}. 

An illustration of the transmitters together with their connectable receivers is shown in Figure [l] 
The family y W can be used to express a variety of frustration events for the transmitter Xi. 
For instance {yW = 0j describes the frustration event that the transmitter Xi is isolated, in 
the sense that it fails to communicate with any of the receivers. Similarly, if B r (Xi ) denotes 
the open Euclidean ball with radius r centered at X % , then Y^> C B r (X t ) encodes the event 
that Xi can only communicate with receivers at distance at most r. 

Before we state our first main result, let us introduce the precise assumptions on the trans¬ 
mission powers and fading variables. We assume that the transmission powers {Px} x eWL d form 
an iid random field whose existence is guaranteed by Kolmogorov’s extension theorem. Note 
that only the subset of powers {Pxf\i> l is relevant, but it is notationally convenient to work 
with the random field indexed by the full space M rf . A similar remark holds for the random 
fading field {F x ,y} x ,y£R d - It can reproduce two different kinds of fading effects. First, a contri¬ 
bution stemming from a suitable random environment such as slow fading, which is typically 
spatially correlated. Second, effects such as fast fading, that are idiosyncratic to the pair (x, y) 
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and therefore do not exhibit spatial correlation. To be more precise for the first contribution, 
we assume Z to be a homogeneous Poisson point process with intensity Ae > 0 modeling the 
random environment. Moreover, we use an iid random held {U x ^ y } xy£ ^d consisting of random 
variables uniformly distributed on [0,1] for the idiosyncratic effects. Then the random fading 
held can have the following general form 

— T(i/ X, Z X , U X) y) 

where $ is measurable and positive. In particular, the construction is such that the fading 
held is spatially translation invariant, i.e., {F x+ZjV+z } is equal in distribution to {F x , y } for any 
z <E R d . 

The dependence of <J> on its second component should be local in the sense that there exists an 
increasing function s env : [0,oo) —> [0,oo) such that $>(z,tp,u) = n I? Senv (i 2 i)(o), u), where 

-^senv(|z|)(°) denotes the Euclidean ball of radius s en v(|z|) centered at the origin. Moreover, 
letting U be a single uniformly distributed random variable on [0,1], we assume that there exist 
N > 0, s max > s m i n > 0 such that for any z € and any locally hnite ip C the distribution 
function q z ^ : t (->• P(<f>(z, ip, U)^ 1 < t) 

(i) is globally Lipschitz with Lipschitz constant N, 

(ii) q x>ip (s) = 0 for s < s min and q x , v (s) = 1 for s > s max . 

The second condition ensures that the fading variables have support bounded away from zero 
and infinity. We assume the same for the power variables P x . Moreover the random objects X, 
Y, Z, {P x } and {U x>y } are independent. 

We provide an example illustrating possible fading helds within the above framework. For in¬ 
stance a Boolean model S = (J z ez Bi(Zi) can be interpreted as randomly distributed obstacles 
in a city. If the line of sight between transmitter x and receiver y is blocked by some building 
the signal propagation is diminished. That is, F X)V = exp ( — 1 [ x ,y] r E.^<n) J 1 (U Xty ) where J is a 
globally Lipschitz distribution function of a random variable which is bounded away from zero 
and infinity. We note that for modeling urban environments it is important to take into account 
the effects of correlated fading variable due to fixed obstacles. See also 0. 

Our first main result will provide a large-deviation principle (LDP) for the empirical measure 
of the family of all connectable receivers — X t such that Xy is contained in the box A n = 
[—n/2, n/2] for large n. To make this precise, we first note that each Y ® is a random variable in 
the measurable space (Nf,A/f). Here Nf is the family of all finite subsets of that is endowed 
with the (j-algebra A/f generated by maps of the form ev# : <p i->- 0 B), for any Borel set 

B C R d . In fact, N f is also a Polish space, see [9, Section A.2.5]. Now, knowing the distribution 
of the empirical measure 

Ln = Ta T S Y(i)-Xi 

1 n 1 Xi£A n 

we can answer questions such as: 

• What is the probability that, when spatially averaged, a certain proportion of trans¬ 
mitters in A n are isolated? 

• What is the probability that, when spatially averaged, a certain proportion of trans¬ 
mitters in A n have l receiver in an r proximity? 

Apart from these examples, L n can be used to describe more general events like an average 
number of connectable receivers per transmitter, i.e., |A n | _1 X^x-eA n • 

The empirical measure L n is a random variable with values in the measurable space 
(Adf(Nf), H cy (Adf)). Here, Adf(Nf) denotes the family of all finite measures on Nf and B cy (Mf) 
is the er-algebra generated by the evaluation maps p > p(B), where B is any bounded Borel set 
of Nf. Since our first main result provides a level-2 LDP, the r-topology on Aif will play an im¬ 
portant role. This topology is generated by the maps p (-)• p(B) where B is any bounded Borel 
set of Nf. We refer the reader to [TQJ Section 6.2] for a detailed discussion of this topological 
space. 
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The LDP allows us to quantify the decay of probability for events away from their ergodic 
limit on an exponential scale. The exponential rate of decay to zero is proportional to the 
volume and the proportionality factor is called the rate function. In order to identify the LDP 
rate function, we first recall the notion of specific entropy of point marked random fields. We 
follow the presentation in [16] and refer the reader also to [ 151 Chapter 15] for further details. 
Let P be a Polish space and write £ for the corresponding Borel u-algebra. Furthermore, let 
N e denote the family of all configurations p C x E whose projection to is injective and 
with image forming a locally finite set. The space N^i is endowed with the smallest er-algebra 
for which all evaluation maps p i— > ff(pC\ (B x F )) are measurable for any Borel sets B , F of 
and E, respectively. Any probability measure on (N^A/#) is called Li-marked point random 
field. Let n > 1 and P be an Li-marked point random field whose realizations are contained 
in A n with probability 1. Moreover, let Q be another Ll-marked point random field that is 
absolutely continuous with respect to P, where / denotes the respective density. Then, the 
specific entropy H(Q\P ) of Q with respect to P is defined as 

H(Q\P) := P(f log /), 

where P(/ log /) denotes the expectation of / log / with respect to P. This definition is extended 
to random point fields Q that are not absolutely continuous with respect to P by putting 
H(Q\P) = oo. Finally, if P, Q are any P-marked point random fields we introduce the notation 

h(Q\P) :=sup-’-Lf(Q A jP An ), 

n> 1 | A n| 

where P Ari , Q\ n denote the projection of P, Q to A n . 

In the following, we write Vg for the family of all stationary P-marked point random fields 
of finite intensity. Here, a stationary P-marked point random field is a probability measure on 
Ng that is invariant with respect to shifts on M. d . The intensity of Q is defined as 



#{(xi, a) G p : Xi G [0, \} d }Q(dp). 


We also need the notion of the Palm version of a stationary point random field as defined for 
example in [2j3J- The (unnornralized) Palm mark measure Q° associated with Q G Vg is given 
by 

Q°(F) = [ #{(xi, e t )Gp: (x i} <*) G [0, l] d X P}Q(d<^), P G £. 

JN E 

In other words, after normalization, Q° describes the distribution of the marks of Q. 

The concept of random marked point random fields is very flexible so that the probabil¬ 
ity space associated with X, Y, Z, {P x } and {U X:V } can be encoded in this framework, see 
Section [2721 for details. 


Let us state the first main result of this paper, an LDP for the empirical measure of con¬ 
nectable receivers associated with transmitters in a large box. Starting from a stationary point 
random field Q of transmitters, receivers and environment, we define Q* as the Palm mark 
measure of the stationary Nf-marked point random field defined by {(Xj, — LQ)}j>i. 


Theorem 1. The random measures {L n } n > i satisfy an LDP in the r-topology with rate |A n 
and good rate function 

T(Q) = inf h(Q |P). 

QeP e 
Q *=Q 

That is for all A G B cy (JA f) 

lirnsup —!— logP(L n G A) < — inf _T(Q) 
n—►oo |A n | QgA 


and 


lim inf ——-logP(P n G A) > — inf 1(Q) 
ra-> oo |A„| Q&A° 
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where A denotes the closure and A° the interior of A respectively. Moreover, the function X is 
lower semi-continuous and has compact level sets. 

To prove Theorem [TJ we make use of the level-3 LDPs established in m (see also mm 
for related results). However, the long-range dependencies induced by the interferences prevent 
us from applying the contraction principle directly. Similarly to [T|, we first have to perform a 
truncation step and consider an approximate model with finite-range dependencies. In order to 
deduce Theorem [l] from the level-3 LDP in the truncated scenario, we show that by a suitable 
choice of the truncation range, the truncation error becomes arbitrarily small. 

In our second main result, we investigate how the connectable receivers associated with a 
typical transmitter located at the origin behave as the connection threshold t tends to zero. 
Since this scenario turns out to be more complicated than the one considered in Theorem [lj we 
impose stronger additional assumptions. To be more precise, we assume that £{r) = min{l, r~ a } 
for some a > d, the transmission power at the origin is fixed (say equal to 1) and that there 
is no random environment Z. That is, {F x ,y} x ,ye R d are iid aR d we put <?( a ) = ^{F~y < a). 
Moreover, we assume that there exist N > 0 and s m j n > 0 such that 

(i) q is globally Lipschitz and globally Lipschitz in its first derivative, both with Lipschitz 
constant N, 

(ii) q(s) = 0 for s < s m j n and q(s) > 0 for s > s m j n . 

To begin with, we provide some important preliminary observations: First, we note that the 
receivers connectable to the origin, namely 

Y l = {Yj € Y : SINR(o, Yj,X U {o}) > t}, 
form a Cox point process with random intensity measure Mf given by 

M t (B) = A r [ T(t- 1 £(\y\),y)dy, 

J B 

where 

r(a, y) = E(q(aI(y)~ 1 )\X) for a > 0 and y e R d . (1) 

In other words, T is an expectation with respect to the fading field in the interference. More 
precisely, it is the conditional expectation on the transmitter process X = {(Xj,Pj)} carrying 
also the transmission powers as marks. For instance, this observation implies that the prob¬ 
ability for the origin to be isolated is given by pt = E exp(—Aft(M d )) and tends to zero as t 
tends to zero. The representation of the isolation probability provides a strong hint that the 
Varadhan-Laplace technique from the theory of large deviations (see e.g. m) could be a use¬ 
ful tool in the analysis of the asymptotic behavior of pt as t tends to zero. In particular, pt 
should decay exponentially as t tends to zero. The exact form of this decay is presented in 
Corollary [3j In Theorem [2j we give a more general result describing the exponential decay of 
unlikely numbers of connectable receivers in space. Throughout the entire manuscript, /3 = 1/a 
denotes the inverse of the path-loss exponent. Furthermore, we put A[ = A 2 r WSm . mt ^-p, so that 
g(t _1 SINR(o, y)) = 0 if y 0 A' t . In the following, we write Pois s for the stationary point random 
field induced by a homogeneous Poisson point process with intensity s > 0. 


Theorem 2. The random measures {|A(| 1 Y t (t l satisfy an LDP in the weak topology 

with rate |A(| and good rate function given by 


where 


%) = 


Im ( ip(y))dy if dip/dx = ip exists, 


oo 


otherwise, 


Zy( s ) = + /i(Pois a |Pois ARQ(r( |j,|-a o)) )), (2) 

and <Q)(r(|y| -Q , o)) denotes the expectation of T{\y\~ a ,o) for a stationary marked point process 
X = {(Xi,Pi)} that is distributed according to Q. 
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Note that in contrast to Theorem [I] the probability measures Q G Vg in ([2]) are distributions 
only of the transmitters X and their transmission powers P. Setting ip = 0 gives the decay of 
isolation probability, this is the content of the following corollary. 

Corollary 3. 

= lim|A^| _1 logEexp(- Ar f T(t~ 1 £(\y\),y)dy) 

= - [ inf (h(Q|P) + A R Q(r(|y|-“,o)))dy. 

J A' 

Large-deviation principles in SIR-based networks have already been considered in mm- 
However, the question treated in Theorem [2] is in a certain sense dual to the ones discussed 
in [T2I [23] . In those papers a large-deviation principle was derived for the interference at the 
origin caused by the signals from other users. We investigate a scenario where the origin sends 
out a signal and we are interested in the interference at the location of the other users. 

The idea for the proof of Theorem [2] is to introduce a stationary point process that carries 
more information than Y t . For this point process, we first establish a level-1 LDP based on the 
results of HE and then deduce a path-space LDP using the Dawson-Gartner technique. The 
proof is concluded by an application of the contraction principle. 

Corollary [3] shows that pt decays exponentially in t~a and provides a variational characteri¬ 
zation of the rate function. However, for the purpose of estimating the actual value of pt , our 
asymptotic result has two drawbacks. First, in Corollary [3j we do not make any claims as re¬ 
gards to how small t should be for the asymptotic to be an acceptable approximation. It is not 
at all clear from the variational formula how to compute (or even approximate) the asymptotic 
rate function. Nevertheless, when estimating the isolation probability pt via Monte Carlo sim¬ 
ulations, our large-deviation result can be used to devise an importance-sampling scheme that 
substantially reduces the estimation variance. In the field of stochastic processes, large-deviation 
techniques have emerged as a powerful tool to find suitable importance-sampling densities [21 
Chapter 6.6], but so far have not found widespread use for spatial rare-event problems. 

As a notable exception, we mention [22], which deals with rare events arising from large values 
of the interference measured at the origin. In that paper, it is shown that the asymptotically 
efficient importance-sampling density is given by a certain inhomogeneous Poisson point process. 
In our setting, the variational characterization in Theorem [2] suggests that the asymptotically 
optimal density is not given by a Poisson point process, but by a collection of location-dependent 
Gibbs processes. Still, in a first step, we provide simulation results illustrating that using an 
isotropic Poisson point process already leads to substantial variance reduction. Let us also note 
that importance sampling for Gibbs processes on the lattice has been studied in J6|. 

The present paper is organized as follows. In Sections [2] and [3] we provide the proofs for 
Theorems [l] and [2] respectively. Section [3] also contains the proof of Corollary [3] Finally, in 
Section [4] we describe two importance-sampling schemes and provide some simulation results. 


2. Proof of Theorem Q] 


As mentioned in Section [lj in order to prove Theorem [lj we use the classical level-3 large- 
deviation result for Poisson point processes [151 . Theorem 3.1]. However, the interferences induce 
long-range interactions that are not immediately compatible with the topology T£ of local 
convergence that is used in [ IB] . To resolve this issue, we will proceed similarly to [I] and show 
that a suitable truncation of the path-loss functions appearing in the interference expression 


induces only a negligible error, see Section 2.1 After this truncation, we show in Section [2.2| how 
the LDP for the stationary empirical field [161 Theorem 3.1] can be used to prove Theorem [I] 


2.1. Truncation of the path-loss function. First, we show that only an asymptotically 
negligible error occurs when disregarding transmitters close to the boundary of A n . This is 
a well-known consequence of the Poisson concentration property [8] Chapter 2.2], but for the 
convenience of the reader, we provide a detailed proof. 
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Lemma 4. Let b,£ > 0 be arbitrary. Then, 

1 


lim , . 

n—t oo 


logP(X(A n \ A n - b ) > e|A n |) = -oo. 


Proof. Let <5 = Aj(l — (1 — b/n ) d ), rn = 5n d and r = en d , then the Poisson concentration 
inequality [8j, Chapter 2.2] implies that 

P(X(A n \ A n _fc)) > r) < ( m/T) T e r ~ n = {6e~ 1 ) end < exp(en d log^e' 1 )). 

Since log(e<5e _1 ) tends to —oo as n —> oo, this proves the claim. □ 


Next, we show that truncating the path-loss function in the interference at a finite threshold 
only leads to a small error provided that the threshold is chosen sufficiently large. To be more 
precise, for 6 > 1 we put £b( r ) = ^(0 if r < b and £b( r ) = 0 if r > b. Furthermore, we define 


I b {Xi,y) = w + Y, p x 3 F X j, y e b (\Xj - y\), and S\m b {Xj, y) 


PXjFx^edXj-y]) 
h(Xi,y) 


and 


Lt 



^Yih.b-Xi- 

Xi^An 


where Y^ ,b = {Yj £ Y : SINR;,(Xj, Yj) > t} denotes the point process of b-connectable receivers 
for the transmitter Xj. We show that when using the total variation distance 


d JS/ (L n ,L b n ) = sup \L n (B) - L b n (B)\, 

B&Xf 


the random measures } n>1 are exponentially good approximations of the random measures 
{Ln} n >i in the sense of jlOl , Definition 4.2.14]. 


Lemma 5. Let e > 0 be arbitrary. Then, 


lim lim sup 

b^-oo n—>-oo 



logP(d TV (L n ,L^) >s) 


—oo. 


Proof. To be specific and for notational convenience let us assume that the support of the power 
variables is contained in [s~a X , * ]. The definition of the total variation distance implies that 

d T v(L n ,L b n ) < e An = Y® + Y^ b }. 

Next, by Lemma [4j we only need to consider those A* that are contained in A„_ 2 ro , where 
ro > 0 is chosen such that £(ro) < wsY m t. Then, almost surely, for Xj £ A n _ 2 ro and Yj £ A£, 
SINR;,(Xj, 1}) < t for all b > 1. Consequently, it suffices to bound the number of transmitter- 
receiver pairs (. Xj,Yj ) £ X x Y such that Xj £ A n _ 2 ro , Yj £ A n and SINR(Xj,Y)) < t < 
SINRb(Aj, Yj). In fact, it suffices to focus on the receivers in these pairs. Indeed, let us call Yj 
b-pivotal if there exists some transmitter Xj such that the pair (X l . Yj) has these properties. 
Then, since we assumed that q x ,tp(r) = 0 for r < s m j n , for each receiver Yj there exist K = 
|"t _1 s^ ax /sA n ] transmitters A(Yj,X) = {X; 15 ..., Xj K } such that SINR^X.;, Yj) < t if Xj 0 
A(Yj,X). Hence, it suffices to show that for every e > 0, 

lim lim sup —r —- logP(#{Y) £ A n : Yj is 6-pivotal} > e|A n |) = —oo. 
b^oo n—Hx> |A n | 

In order to do so, we use the exponential Markov inequality with s > 1 and estimate 
P (#{Yj £ A n : Yj is 6-pivotal} > e|A n |) 

< exp(—se|A ra |)Eexp(.s^{y ; ' £ A n ■ Yj is 6-pivotal}). 

Hence, it suffices to show that for every s > 1, 

lim lim sup ——- logE exp(s#{Y) £ A n : Yj is 6-pivotal}) = 0. 
b^foo n —>oo |A n | 


7 









The point process of receivers that are 6-pivotal form a stationary Cox point process with 
random intensity measure 


M'(B) = Ar / P(y is 6-pivotal|AT, Z)dy, B G B(R d ), 

J B 

where we think of X = {(A*, Pi)}i> 1 as a marked point process and we evaluate the probability 
with respect to the fading variables associated with the pairs (y, Xj)i>i. Since q is assumed to 
be globally Lipschitz with constant N, we arrive at 

P(y is 6-pivotal|X, Z) 

< nS\NR(X h y) <t<S\m b (Xi,y)\X,Z) 

Xi&A(y,X) 

< J2 e [PxA\Xi-vMXi,vT\ PxA\Xi-y\)h{X u y)- l ]\X,Z) 

Xi£A(y,X) 

< AXi-yWsAj-'w-^AXi^-AXi^iXiZ), 

Xi£A(y,X) 

which is at most 

SJ2AXi-y\)-l b (\Xi-y\) 

i> 1 

where 5 = K£(0)Ns m ? n t~ 1 w~ 2 . In particular, we obtain that 

M'{B) < S' [ Z(\ x i - y\) - h(\Xi - y|)dy, 

Jb i> i 

where 5' = Ar5. Hence, using the formula for the Laplace functional of a Cox point process, 
we get that 


Eexp [s#{Yj € A 2 n : T, is 6-pivotal}] 

< Eexp [(e s - 1)5' [ ^ " V\) ~ Xi ~ y\)M 

■1^271 j>! 

= exp [A t f exp ((e s - 1)5' [ £(\x - y\) - £ b (\x - y\)dy) - ldx]. 
</K d J A2 n 

Notice, that we can bound the integral 

/ exp [t £(\x - y\) - £ b (\x - y\)dy] - ldx 
JR d JAor, 


' A 2 n 

< / exp [t 


■£(|a; — y|)dy] — ldx + / exp [r / £(|x — y|)dy] — ldx. 

jR d \A 4 „ J A2n 


where r := ( e s — 1)5'. In the next step, we derive bounds for these expressions separately. For 
the first, we get that 


f exp [r f £(\x - y\)dy\ - ldx = 4 d [exp (r f £(\y\)dy) - l] 

J A 4 „ JR d \B b (x) JR d \B h (o) 


1 

'A 4n 


R d \B b (o) 


which tends to zero as b tends to infinity. For the second expression, we note that for x € 
and y G A 2n , 


lx _ \ = \z~y\ + XAx~v\ > I*-i/I + \v\ > (1 

1 1 + Vd “ 1 + \fd ~ 


x . 



Consequently, using that £(r) E o(r d ), we have for large n 

1 f 


nr 


< — 
n‘ 


< 


exp [t / £{\x — y|)dy] — ldx 

* Mn 

—7 [ exp \T(2n) d £((l + \/d) _1 |x|)l — ldx 

jRd\ 


' E d \A 4 


[ r2 d+1 £({l + Vd)~ 1 \x\)dx = 

J E d \A,i n 

which tends to zero as n tends to infinity. 


r2 rf+1 (l + Vd) d £(\x\)dx, 


4n(l + -\/d) - 


□ 


2.2. Application of LDP for the stationary empirical field. In order to apply [16, The¬ 
orem 3.1], we need to relate the empirical measure of connectable receivers to the stationary 
empirical field considered in [16]. Here the first task consists in encoding the probability space 
carrying the point processes of transmitters X, the point process of receivers Y . the random 
environment Z, the transmission powers {P x } and the iid family {U XtV } in the framework of 
stationary marked point processes. To be more precise, we put X = {E, R,T} and consider 
the mark space E = X x (0, oo) x [0,1] N equipped with some complete and separable metric. 
Furthermore, we let V denote an independently E-marked homogeneous Poisson point pro¬ 
cess with intensity X a . The mark distribution on E is a product of three distributions 

defined on the spaces X, (0, oo) and [0,1] N , respectively. First, on X, we choose the distribu¬ 
tion which assigns a E X the probability A (r /(^ (T / gS A CT /). Second, on (0, oo) we choose the 
distribution of the transmission power P x considered in Section [TJ Third, the distribution on 
[0,1] N describes a family of iid random variables that are uniformly distributed [0,1]. The 
Poisson point process Z that generates the random environment is represented by elements of 
V = ( Vi,cri,Pi , with cjj = E. Elements of V = ( Vi,Oi,Pi , (C^,j)j>i)i>i with cq = T 

are thought of as transmitters and are denoted by X. Elements of V = (Vi,<Ji,Pi , {Ui,j)j>\)i>\ 
with <7j = R are thought of as receivers and are denoted by Y. We note that the power variables 
Pi have not meaning if cq ^ T. The random variables U XjV should be thought of as being 
attached to the the transmitters. Moreover, proceeding as in ESI Section 1], let 


P per ’ n = U ((FnA„) + n S ) 


denote the periodic spatial continuation of V n A n . The stationary empirical field is defined as 


Rn,v '■= 


1 


|A. 


1 


ypei 


i.- v dv 


n\ J A„ 


where !/ per ’ n — v = {(vj — v,ej)}j>\ is the spatial translation of l/ per,n by v. Now, we let 
yper,n,6,(i) denote the family of periodized receivers that have a 6-connection to the transmitter 
Xf er,n = (xi, T, Pi, More precisely, 


yper,n,b,(i) _ | y. 




where 

P X i,yj ^ (l-Jj x ij Z P •^iiPi,'\’(yj—Xi,YP e, ’ n — X i))i 

and where the integer T (yj — Xi , Y per,n — Xi ) > 1 is defined as follows. If k > 1 is such that yj — Xi 
is the /c-th closest element in y per ’ n — x% to the origin, then we put d>(yj — Xi,Y per,n — Xi) = k. 
This construction will ensure translation invariance for the periodized version. The empirical 
measure L^ r ' b of 6-connectable receivers associated with transmitters in A n when the network 
is based on periodized configurations can also be expressed as a function of R n ,v- Indeed, by 
the same technique that was used to define the individual empirical field in [16], we arrive at 



I A, 


X 


A^xper,n,b,(i)_^ 


xieXP er > n nA„ 




</(P per ’ n 


u)du, 
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where 

g'iy P«r.n_„) = £ 5 (yp-,n_ x . )i 

Xi-vG(X P er - n -u)nAi 

and where g is the Dirac measure concentrated on the family of 6-connectable receivers from the 
origin multiplied with the indicator function that the origin is a transmitter. Next, we prove 
that the random measures {L^ er,fe } n >i and {L b n } n >i are exponentially equivalent (in the sense 
of |TOs Definition 4.2.10]), when using the total variation metric. 

Lemma 6. Let e > 0 be arbitrary. Then, 

lim — 1 — logP(d T v(^ er,6 ,-bn) > <0 = 

n—>oo |A n | 

Proof. As in Lemma [HJ choose ro > 1 such that £(\x — y\) < ws^- in t if \x — y\ > ro- In 
particular Y^ ,b C B ro (Xi). Further by the truncation of the interference, to decide whether 
Yj E yW’ 6 ^ suffices to look at transmitters in BbiYj). As a consequence, the family Y^ ,b 
depends only on the network configuration in B ro+b+Sem ^(Xi) . Hence, djv(L^ er ’ fe , L b n ) < ff(Xf\ 
A n \ A n _ 2 ( r - 0 _ ) _fe_|- Sen v(b))), and the claim follows from Lemma |4j □ 

Now, we are in a position to provide a proof for the LDP asserted in Theorem [l] when 
L n is replaced by . Let Q be the distribution of some stationary .E-marked point process 
V = Then, we define Q* ,b as the Palm mark measure of the marked 

point process (Xi,Y^ ,b — Xf). Here, as above, X, E V are interpreted as transmitters and 
y(d> 6 c V as the 6-connectable receivers. 

Proposition 7. The random measures {L^} n >i satisfy an LDP in the r-topology with rate |A n | 
and good rate function 

l b :Q^ inf h(Q IP). 

QePe 

Q*’ 6 =Q 

Proof. First, we note that the map <!>{,: Q >—)■ Q* ,b is continuous with respect to the r-topology. 
Indeed, for Q n —> Q and B E A/"/, the locality that is established after truncating the interfer¬ 
ences gives that |Q* ,fc (E) — Qn’ 6 (E)| —> 0 as n —> oo. As &b(Rn,v) = L% er ' b , we can apply [Tbl 
Corollary 3.2] and the contraction principle. Thus the random measures {L£ er,b } n >i satisfy an 
LDP with good rate function lb. Finally, Lemma i shows that {L£ er,b } n >i and {L b n } n >i are 
exponentially equivalent with respect to the total variant distance. This implies exponential 
equivalence of {Lj) er,b } n >i and {L b n } n >\ when evaluated on an arbitrary Borel subset of Nf. So 
the claim follows from m Corollary 1.10, Remark 1.4]. □ 

The same arguments also prove the following result, where we consider the marked point 
process (JQ, Y^ ,b — Xi, Y^ ,b ' — Xf) at different truncation thresholds b' > b > 1. Starting from 
Q E Vg, the associated Palm mark distribution is denoted by Q* ,b,b . 

Lemma 8. Let b' > b > 1. Then, the random variables { | A ’ <fj z {X t E A n : Y^X / Y ^’ b }} n >i 
satisfy an LDP with rate \ A n \ and good rate function 

s i —y inf /i(Q|P). 

Finally, we complete the proof of Theorem [I] In Lemma [HJ we showed that {L^} n >i are 
exponentially good approximations of {L n } n >i and hence an application of eh Theorem 1.13] 
is natural. 

Proof of Theorem [7} As mentioned in the previous paragraph, eh Theorem 1.13] implies that 
it suffices to verify the following condition. For every e, K > 0 there exists b > 1 such that 

sup d T \/(Q*,Q*’ b ) <e. 

Q&Vg 

h(Q\F)<K 
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We show that a slightly stronger statement holds, where djvCQ*: Q*’ 6 ) is replaced by <Q)*(y(°)> 6 ^ 

y(°>). 

Lemma [5] shows that there exists bo > 1 such that if b' > b > bo, then 

limsup|A n |- 1 logP(#{W G A n : Y^ b ' ± Y^ b } > e|A n |) < —K. 

n—>00 

Hence, the LDP from Lemma [8] yields that 

inf h(Q IP) > K. 

Qo &Vg 

<Q*’ b ’ b ' (Y(.o), b ' b )>e 

In particular, if /i(Q|P) < K, then Q*’ &, A(y(°)A' ^ y(°)’ b ) < g. as required. □ 

3. Proof of Theorem [2] 

The difficulty in proving Theorem [2] is that the connectable receivers associated with the 
origin are not stationary, so that we cannot use LDPs for the stationary empirical held directly. 
Therefore, we first consider a more general stationary marked point process from which the 
connectable receivers can be reproduced by an application of the contraction principle. Since we 
need a path-space LDP for this stationary marked point process, we proceed as in the classical 
proof of Mogulskii’s Theorem m Theorem 5.3.1] and use the Dawson-Gartner technique to 
deduce the path-space LDP from the finite-dimensional marginals. 

In order to define a suitable auxiliary stationary marked point process, we consider the 
random measure 

r roo 

M\-) = A r / / l{-}t'j / (ds)dy, 

J A' Jo 

where v y ([0, s]) = T(s,y) = E(( 7 (s/(y) _1 )|X) see also ([!]). Then, we let Z l = {(Ay, *Sy)} denote 
a Cox process with this random intensity measure and define the two-parameter held Y* ,t = 
{y ’ {x , s)}(a; 1 s)£A , 1 X [0,00) by 

Y*’\x,s) = Z t { A t (s) x (0, s]), 

where 

d 

i=i 

In particular, for any hxed (x, s) G A' x x [0, 00 ), conditioned on X the random variable Y* ,t (x , s) 
is Poisson-distributed with parameter Ar f At ^T(s,y)dy. Moreover, Y*J is a random variable 
with values in L- mc , the space of [0, oo)-valued, bounded and coordinate-wise increasing functions 
on A] x [0, 00 ). 

In the following, we put /xq(s) = ArQ(P(s,o)) and note that the derivative exists 

since q is differentiable and ^q(s) is Lipschitz continuous with Lipschitz constant N. 

Similar to |101 Section 5.3] we introduce the notion of absolute continuity for increasing 
functions F : A] x [0, 00 ) —>• [0, 00 ), (x, s) <—> F(x,s). For the convenience of the reader, we 
reproduce some of these definitions and observations. F defines a additive set-function on the set 
of cubes. More precisely for A = (ai, b\] x ( 02 , 62 ] x • • • x (ad, bd] x (a^+i, &d+i] we will sometimes 
write F( A) := Ylu a ( u )F(u) with a(u) := (— l) p where p = #{k : Uk = afc} and the summation 
extends of all corners a of A see m Chapter 3]. It follows from Caratheodory’s extension 
theorem, that any right-continuous F G L\ nc induces a unique measure fij?on A\ x [0,oo) with 
the Borel sigma-algebra satisfying 

d 

MllH^Af] x [ 0 , s ]) = F(e,...,t d ,s), 

2—1 

for any s > 0 and (^,..., £ d ) G A',, see Theorem 3.25]. 
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The function F is called absolutely continuous if F is right-continuous and fj,p is absolutely 
continuous with respect to the Lebesgue measure on A' x x [0,oo). We write dF/(dxds ) for its 
Radon-Nikodym derivative. Let 

AC'l := {F : F is absolutely continuous, F(x, 0) = 0 and F(— |A / 1 | 1 / d /2, f 2 ,... ,£ d ,s) = 

= F(e, -\K\ 1/d /2, ...,t d ,s) = --- = F(e, • • • ,Z d -\-\K\ l,d /^ s) = 0}. 

In Section 3.3, we will derive Theorem [2] by the contraction principle from the following result. 

Proposition 9. The random fields {|Aj| _1 y*’*(-, satisfy an LDP in the topology of point- 

wise convergence with rate |Aj| and good rate function given by 


1(F) = 


/.vydsfe.'))*/ ifFeACi, 

otherwise, 


oo 


where 

roo 

^*( 5 ) = inf h(Q|P) + / h(g(s)\£n Q (s))ds. (3) 

*£€Pe Jo 

3.1. Finite-dimensional result. In order to apply the Dawson-Gartner Theorem [ TOl Theo¬ 
rem 4.6.1], we first derive the finite-dimensional LDPs. 

Proposition 10. Let —|A / 1 | 1//rf /2 = £0 < £1 < • • ■ < fk < |A , 1 | 1 / rf /2 and 0 = so < s\ < ■ ■ ■ < s r . 
Furthermore, put E = {£o>£i> ■ • • >£fc} and S = {so> si, ■ ■ •, s r }. Then, the random variables 
{(|AJ| _1 F*’*(x, s)), S )gcd x s} t<1 satisfy an LDP with rate |Aj| and good rate function 


ls,s(F) 


1 

h (^\ P ) + £ h (\XFM\ F{AE{x) X S *DI a Mq(' s *)) 5 

r-a’P.d- i=l 1 


where A// Q (sj) = nq(si) - i), and A s (Cu, ■ • • = 1 I'; |(G r G, • 1 ] with (& 1; ... ,& d ) G 

ztd 


The basic idea of proof for Proposition [lO] is to apply the LDP for the stationary empirical 
field [16, Theorem 3.1]. However, in order to cast our problem into a suitable framework, we 
first have to perform a truncation and a periodization step. 

3.1.1. Truncation of the path-loss function. In a first step, we show that truncation of the path- 
loss function gives an exponentially good approximation. Let b > 1, s' > s > 0 and x, x' G Ai be 
such that all coordinates of x' — x are positive. Then, we let Y* ,b,t (x , x', s, s') denote a random 
variable that conditioned on the independently marked Poisson particle process X is Poisson 
distributed with parameter Ar f At(x T b (s', y) — T 6 (s, y ), where T fe (s, y) = E(sL b (y)~ 1 \X), 

d 

A t (x,x') =t~ 13 A(x,x') = t~ 13 Y\( 7T k(x),TT k (x')\, 

i=l 

and 7 Tfc : M. d —> M denotes the projection onto the fcth coordinate. From now on let us again 
assume that the support of the power variables is contained in [0, s m ? n ]. 

Lemma 11. Let b > 1, s' > s > 0, and x,x' G A ; ,. Then, {Y* ,b,t (x, x', s, s / )}b>i,i<i are 
exponentially good approximations of {Y* ,t (A t (x,x') x (s, s r ])}t<i - 

Proof. Conditioned on X , the random variable \Y*' b,t (x, x' , s, s') — Y* ,t (x , x', s, s')| is stochasti¬ 
cally dominated by a Poisson distributed random variable with parameter 

H b = Ar [ T b (s', y) - r(s', y) + T b (s , y) - T(s, y)dy. 

J A t ( x,x') 
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Hence, using the Laplace transform of Poisson random variables, for any a > 1 the exponential 
moment of a\Y* ,b,t (x, x', s , s') — Y* ,t (x, x', s , s')| are bounded from above by 

Eexp(a|y*’ 6,t (x, x' , s, s') — Y*’ t (x, x’ , s, s')|) < Eexp((e“ — 1 )Hf). 

Now, similar to the proof of Lemma [5j H & can be bounded from above by 

A R iV(s + s')s m ^ n w~ 2 [ VVdXi - y\) ~£ b {\Xi - y\)dy, 

J At(x) g 


so that 


Eexp(a|y*’ 6 ’*(.T, x', s, s') — Y* ,t (A t (x,x') x (s,s'])|) 

<Eexp((e a -1)A rN(s + s')s~f n w~ 2 I 'Y'H\X i -y\)-t b {\X i 

V JA t (x,x’) 

Now, we conclude as in Lemma [5] 


□ 


3.1.2. Periodization of the integration domain. Next, we show that replacing the quantity 
Y*’ b,t (x,x',s,s') by a periodized variant is exponentially equivalent. To be more precise, let 
^>l)'S / >s>0 and x, x' £ be such that all coordinates of x’ — x are positive. Then, X per,t 
denotes the periodization of An A t (x,x'), i.e., 


* per>t = U (\^t{x,x')\ 1/d z + XnAt{x,x')). 

z£Z d 


As in Lemma 11 we let Y* ,per,b,t (x, x', s, s') denote a random variable that conditioned on X is 
Poisson distributed with parameter 


Ar [ T per ' b (s',y) -T per ' b (s,y)dy. 

J A t{x,x') 

Here, r per,6 (s ) y) = K(q(sL per ' b (y)~ 1 )\X per ' t ) and I per,b (y) is the interference at y in the periodized 
configuration computed using truncated path-loss functions. 


Lemma 12. The random variables {y*’ per,6 A(x, x', s, s / )} t<1 are exponentially equivalent to the 
random variables {Y* ,b,t (x, x 1 , s, s ')} f<1 - 

Proof. Since we consider truncated interferences, we have that I per,b (y) = I b {y) for all y £ 
Af(x,x'), where Af(x,x') denotes the subset of all y £ A t (x,x') such that B^y) C A t (x,x'). 
In particular, | Y*- per ’ b ’ t (x, x', s, s') — Y* ,b,t (x,x', s, s') | is stochastically dominated by a Poisson 
random variable with parameter 2A R |At(x, x')\Af (x, x')\. Now, we can conclude as in Lemma[6] 
by making use of the Poisson concentration property. □ 


3.1.3. Application of LDP for stationary empirical fields. We have seen that truncating the 
interference and considering a periodization does not have an effect on {Y* ,t (x, x', s, s')}* <1 in 
the LDP asymptotics. Now, we derive an LDP after these modifications have been implemented. 
We put fA(s) = Q(T b (s, o)). 


Proposition 13. The random variables {|A£| 1 ffY* ,per ' b,t (x,x', s, s ')} r< i satisfy an LDP with 
rate |A£| and good rate function 


q-X,x' ,S,s' 

H.n 


(a) 


l A (*>*')! M WW + -/^q( s ))) 


(4) 


Let us recall from [IQi Equations 1.2.12 and 1.2.13] that if the random variable considered 
in an LDP is measurable with respect to the Borel cr-algebra on the underlaying topological 
space, then the proof of the upper and lower bound can be done directly for closed and open 
sets, respectively. We use this in the sequel without further mentioning. 
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We prepare the proof of Proposition 13 by a lemma. First, we note that m Theorem 3.1] 
gives the following auxiliary result, where we put 


r per, 6 (s / , y) — r per ’ fc (s, y)dy. 


M aVjt = M aVjt (x,x',s,s') = A R |A t (x, x')| 1 

J A t{x,x') 

Lemma 14. Let F and G be compact and open subsets of [0,oo), respectively. Then, 
lim sup 


and 


—^ i . / tit log IE exp ( - \A t (x, x')\ inf h(a\M av r )) 

t-+0 x')\ aeF 

< - Q inf^/i(Q|IP) + /i(a|/TQ(s') -//q(s)), 

a£F 


liminf ——-— logEexp ( - |A t (x,x')| inf h(a\M ay r )) 

t—>0 J| cl(E.G 

> - inf h(Q\F) +/i(a|/x{j(s') -fiq{s)). 

U&Fb 

a&G 


Proof. In order to apply m Theorem 3.1], we only need to check that the functions 

Q ^ ^nf h( a \n b Q (s') - Mq(s)) 

and 

<Q> ^ mf h( a \PQ( s> ) ~ Mq(s)) 

are lower- and upper-semicontinuous, respectively. First, note that the map Q i-)- p,^(s') — /Tq(s) 
is continuous in the r^-topology, since r b (-,o) only depends on X via X n Now, we 

conclude by observing that a' i—r inf h{a\a!) is lower-semicontinuous as pointwise infinrunr of 
a two-parameter lower-semicontinuous function over a compact set and a' i —> inf ag c h{a\a') is 
upper-semicontinuous as infimum over a family of continuous functions. □ 


Now, we can proceed with the proof of Proposition [l3| 

Proof of Proposition \1S\ The upper bound for compact F is an immediate consequence of 


Lemma 14, since [2TI Lemma 1 . 2 ] implies that 


P(|A£| 1 Y*' per ' b,t (x, x', s, s') 6 F)< Eexp ( — | A i (cc, rc / )| inf hi 

\ a£F 


|A(ai,ai / )| 


\M- 


a v,r 


The proof of the lower bound is more involved. First, we may assume that G is an interval, 
i.e., G = [ 0 , 7 ) or G = (y_, 7 +) for some 7 , 7 -, 7 + > 0. Next, introduce the function f(k ) by 
/(0) = 1 and f{k) = e - 1 /( 12fc )(^/27r k)~ l for k > 1, and put G* = Zfl (|A^|C?). Then, by [5T[ 
Lemma 1.3], 

P(|Aj| _ 1 y ! "’ per ’ & ’ t (x, x', s, s') € G) 

> Eexp ( - ^inf t - log f{k) + h(k\\A t (x,x')\M av j)^ 

> Eexp |Ai|(|A ;|" 1/2 + |A(x, x')\ Jurf M av , t h(k\A t (x, x , )| _ 1 M"J t ))), 

where h(k\A t {x, x 1 )^ 1 M~\) = h(k\A t (x, x')^ 1 is a short notation. Now, we distinguish 

between the cases where G contains 0 and where it does not. We claim that if G = [ 0 , 7 ) and 
e > 0 , then 

fc hT M aVtt h(k\A t (x, x / )| _ 1 M a “| ]t ) < e + mf M av , f fe( \ A ^\ M ~] t ), 

provided that t > 0 is sufficiently small. Once this claim is proven, Lemma [14] completes the 
proof of the lower bound for the case G = [0, 7 ). Let e > 0 be arbitrary. Then, under the event 
M av ,t < e, we deduce that 









On the other hand, if M av j > e, then for every g E G, 

where k(g) > 1 is chosen as the element of Z n (|A' t |G) such that k(g)\A t (x, x')\ 


-l 


mini¬ 


mizes the distance to g\A(x,x')\ 1 . In particular, uniform continuity of h(-) on the interval 
[ 0 , 7 |A(x,x , )| _ 1 e _1 ] implies that 

inf M av gh(k\A t (x, x')^ 1 M~\) < e + inf M av t h{g\A(x, x , )| _ 1 M a “ 1 t ), 
fceG 4 ’ geG ’ 

for all sufficiently small t > 0. Finally, we deal with the case, where G = ( 7 -, 7 +) and observe 
that if M aV; t > e, then we can conclude as before. To be more precise, 

Eexp ( - \A t (x,x')\ inf M av t h(k\A t (x, 

> Eexp f- |A t (x,x')|(-e + inf M aVit /i(g|A(x,x , )| _ 1 M“^))') 

-El{M aVjt < e} exp ( - |A t (x, x')\(-e + inf M aVjt h(g\A(x,x')\^ 1 M~ 1 t ))). 

geG 

Now, for any K > 1 there exists e > 0 such that |A(x, x')\~ 1 M~^ t ) > I\ if M aVjt < e. 

In particular, 

exp ( - \A t (x, x')\{-e + mfM aV)t h{gM~] t ))) < exp ( - \A t (x,x')\(-£ + K)), 

geG 

which completes the proof of the lower bound. 

Since Y*’ per,b,t (x, x', s, s') is stochastically dominated by a Poisson random variable with pa¬ 
rameter ArIA^Xjx')), the random variables {|A / t | _ 1 y* ,per,6 ’ t (®, x', s, s')}<<i are exponentially 
tight. This implies both goodness of the rate function and the full LDP. □ 


Next, using Lemmawe derive an LDP for the finite-dimensional marginals of T* ,per,b, *(-, •). 
In order to state this precisely, it is convenient to introduce some notation. Let — |A / 1 | 1/,rf /2 = 
£0 < £1 < ■■■ < £fc < \A’ l \ l / d /2 and 0 = s 0 < si < ■ ■ ■ < s r . Then, for x = (£i x , - - -, £i d ) 
and s = Si we put x +j h = (£ii+i,..., £j d +i) and s +: s = Sj+i, where we use the conventions 
£fc+i = |A / 1 | 1 / d /2 and s e+1 = 00 . 

Corollary 15. Let —|A / 1 | 1//d /2 = £0 < £1 < • • • < £fc < |A / 1 | 1 / d /2 and 0 = sq < s\ < ■ ■ ■ < S£. 
Furthermore, put S = {£o;£i) ■ ■ • >£fc} & n d S = {so,si,..., s r }. Then, the random vectors 

{(|A*I ^ ’ ’ ( x ’ x +. s ’ S ’ ,s +>S , ))( 3 . iS ) eS d X 5} < 

satisfy an LDP with rate \ A! t \ and good rate function 


1 t<i 


4,s((o. 


zjze^xs) 


Y. 

xGS' 


|A(x,x+ i5 )| jnf^ h(Q\F) + X) /t ( |A(xX i , B )l 




1=1 


Proof. First, we observe that (Y* ,per,b,t (x, x+,s)),j. g cd defines a family of independent random 
vectors, where we put 

y*- per > 6 > t (x,x +iS ) = (y^ per ^(x,x + , s , s , s+ , 5 )) seS . 

Indeed, this is a consequence of the independence property of the Poisson point process X since 
by the definition of the periodization, Y*’ per,b,t (x, x + =) depends on X only via X nAt(i,x + n). 
Hence, if for any fixed x E F. d we can establish an LDP for Y*’ per,b,t (x, x+=) with a cer¬ 
tain good rate function, then m Exercise 4.2.7] allows us to deduce that the collection 
(Y*’ per ’ b,t (x, x+,H)) Tg cd satisfies an LDP with good rate function given by the sum of the indi 
vidual ones. If \S\ = 1, then the LDP for Y*' per,b,t (x, x+s) is precisely the result of Lemma 
and an inspection of its proof shows that it also extends to the case of general finite S. 

In Lemma 


13 
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we have seen that periodization replaces { Y* ,b,t (x , x', s, s / )} T <i by exponentially 
equivalent random variables. Hence, applying [10L Theorem 4.2.13] completes the proof. □ 
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In order to deduce Proposition [TO] from Corollary [T5j we need to undo the truncation approx¬ 
imations. Before we start with the proof of Proposition 10 it is convenient to derive certain 
continuity properties of /Jq and ^/Xq(s) with respect to b and <Q). The technique of proof is 
similar to the one used in the proof of Theorem [l] 

Lemma 16. Let e,K >0 and s > 0 be arbitrary. Then, there exists b > 1 such that if Q € Vg 
satisfies /i(Q|P) < K, then 

and \4zUn(s) - 4-un(s)\ < e. 


(s) fJ,Q(s) < £, U./I.U, dD 

Proof. We first deal with the part of the statement not involving derivatives. Let b' > b > 1 be 

we see that {|A^| _1 / a , r 6 (s, y) — r &/ (s, y)dy} f<1 satisfies 
'unction 


arbitrary. Proceeding as in Lemma 


14 


an LDP with rate |A(| and good rate 


a i —> 


inf 

QeVp 

ij -q(s)-Uq (s)=a 


In particular, the proof of the lemma is completed, once we show the existence of bo > 1 such 
that 

limsup |A (| _1 logPf f |At| _ 1 r & (s,y) - T b \s,y)dy > e) < -K. 

t->o v JA[ 7 

for all b' > b > bo. Note that 

T b (s, y) - T b \s, y) < Nss^ n w~ 2 £ 4'(l y - X t \) - 4(1 y ~ M)- 


i> 1 


(5) 


Hence, using the formula for the Laplace functional of a Poisson point process shows that for 
any a > 0, 

Eexp (a J' T b (s,y) - T b \s,y)dy > e) 

< exp (A T / exp ( aNs~f n w~ 2 [ £(\x - y\) - l b (\x - y\)dy) - Ida;) . 

V J R d J A' 7 

Now, we conclude as in Lemma [5] 

For the part involving the derivatives, note that the derivative of ixq(s) is given by 
ArQ(/(o) _1 ^P(s, o)). Essentially, this means replacing in the above arguments the expression 
q(sl(y ) _1 ) by I(y)~ 1 -^q(sl(y)~ 1 ). This specific form only comes into play in the estimate ([5]) 
which can be replaced by 

< - abKAr 1 )! + - /‘Ar 1 ! 

< (jVsSmin”’ -3 + JV ”’“ 2s mfn) “ A '‘l> “ Wl* ~ A iD)> 


i> 1 


as required. 


□ 


Corollary 17. Let K > 0 and s > 1 be arbitrary. Then, in the rc-topology, 

(i) as b —> oo, the functions Q H > Pq(s) converge to hq(s) uniformly in {Q : /i(Q|P) < K}. 
In particular, Q i —> /mq(s) is continuous on {Q : ft(Q|P) < K}, 

(ii) as b —> oo, the functions Q (->■ ^/Zq(s) converge to j^Pq(s) uniformly in {Q : /i(Q|P) < 

K}. In particular, Q i —> is continuous on {Q : /i(Q|P) < K}, 

(hi) ifb n ->• oo, Qn -t Q and limsup n _ ) . 00 /i(Q n |P) < K, then ^(s) hq(s). 

Proof. Since the first two items are an immediate consequence of Lemma [l6j we only deal with 
the last item. Here, we use the decomposition 

\TQ (s) - /*<£(»)| < Imq(s) - Tq(s)\ + \y b Q (s) - n b Qn (s )| + |/ij£ n (s) - /iQ n (s)|, 
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and conclude as before. 


□ 


Now, we have completed all preparations for the proof of Proposition 10 


Proof of Proposition \10\ In Lemma 11, we have seen that truncation leads to an exponentially 
good approximation. Therefore, combining Corollary 15 with |10, Theorem 4.2.16] shows that 
the random vector { (| A£|~ 1 y* ,t (x, s))^. s ) e =d Xi g} t<1 satisfies a weak LDP with rate function 

-Ek,s(( a x,s)(x,s)) = sup lim inf inf I s,s(i a x,a){x,a))- 

212 >1 0 \ a x,s)(x,s) • I ( a x,s)(x,s) (^ l 'X,s)(x,s) 100^771 

Since the random vectors {(|A£| _1 Y* ,t (a;, s)) ^ s ^ e ^d xS } T<1 are exponentially tight, the proof is 
completed once we show that X4 s(( a x,s)( x ,s)) = 2s,s((ox,a)(x,s))) where 

r 

l~,s((a x ,s)( x ,s)) = I A 0d x +,s)I Q inf g HQ |P) + ^ h ( \A^xi s )\ 1 A Pq(^))- 

x£B. d i =1 

First, we note that q(sI b (o)~ 1 ) is decreasing in b and converges to q(sI(o)~ 1 ) as b —> oo. Hence, 
for any Q € To and x G E d , 

ji™ iz h ( ia(S;. B )i i ■ A ^ (si) ) = s h ( ia(S+, 5 )i i ■ 

2=1 2=1 

In particular, Z4 s(( a x,s)(x,s)) — ,s((flx,s)(x,s))- For the other direction, fix S > 0 and x G 

Then, for each m > 1 choose a sequence (& m ,n)n>l such that lim^^oo b mjn = oo and write 


lim inf 

n—>oo QePfl 

( a 4,s)(x,s) : l( a x,s)s — (®a;,s)s|o °<TTl 1 


.-1 i =1 


= lim inf inf 

6—^00 

( a a;,s)(a;,s) : l( a i,s)s — ( a ^,s)s |oo<2?2 1 


,-i i=l 


Next, for each m, n > 1 choose Q x ,m,n G Pe, and for each m, n > 1 and s G S' choose a' x smn G 


[a X)S — 1 /m,a XjS + 1/m] such that 


inf 

( a J,s)s : | ( a J,s)s — (0'X,s)s\oo<7Tl 1 


: |p) + I> 


2=1 


|A(x,x+, 3 )l 


A^' n (si)) 


> - 


6 + /l(Qx, m , n |P) + X^(|A(x,x+, s )l I A ^«£’,l,n ( S *)) ‘ 

2=1 


If lim sup, 


1 1P) = oo, then 


n—»oo 

lim inf 

n->oc QePe 

( a x,s)(x,s) : l( a a:,s) s — ( a a;.s)s|oo <Til 1 


+ i A ^p ’ ( Si ))—°°) 

2=1 


which is certainly at least as large as Ie,s{(o-x,s)( x ,s))- Otherwise, after passing to a subsequence, 
we may assume that lin^^oo Q x ,m t n = Q x ,m for some Q Xtm G Vq by sequential compactness. 
Furthermore, we may also assume for each 1 < i < r, that lim n ^. 00 a' x mn = a' x m for 
some a! x s . m G [a XjS — 1/m, a X)S + 1/m]. In particular, lower semicontinuity of h implies that 


lim inf ,. 


x,m,n P) > 


/P). Moreover, by Corollary 17 




Hence, another application of lower semicontinuity gives that 


lim inf h 

n—>oo 


x,s^ ,m,n 

|A(x,x+, s )| 


>>) a 
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Arguing as above, we may assume that Q xm converges to some <Q> X as m —> oo. In order to 
conclude the proof of the proposition, it therefore suffices to show that 


lim inf h 

m—>oo 


x.sa ,m 


|A(x,s + ,e)| 




A final application of lower semicontinuity completes the proof. 


□ 


3.2. Application of Dawson-Gartner &; identification of rate function. In Proposi- 
we have shown that the finite-dimensional distributions of the random fields {Y* ,t }t<i 


10 


tion 

satisfy an LDP and we have also identified the good rate function. Hence, the Dawson-Gartner 
Theorem m Theorem 4.6.1] implies that the random fields i satisfy an LDP with 

respect to the topology of pointwise convergence and that the good rate function is given by 

1(F) = supX HiS (F), 

s ,s 


where the supremum is over all finite S C [0,oo) and S C [—|A / 1 | 1 / d /2, |A / 1 | 1 / rf /2]. The proof of 
Proposition [ 9 ] now amounts to showing X(F ) = 1(F). This can be done using an adaptation 
of arguments appearing in the classical derivation of Mogulskii’s Theorem provided in m 
Theorem 5.3.1]. For the convenience of the reader, we provide some details. 


Proof of Proposition^ First assume that F e ACq and let / = dF/(dxds) denote the density 
of F. By non-negativity of Is t s(F), we can assume fk = A / 1 | 1 ' /d /2. Note that X^^s(F) is of the 
form Ylxezd |A=(x)| inf(j f(F(x),Q) where / is convex in the pair (F(x),(Q)) by linearity of pq. 
Hence, also G(F(x)) = iniqf(F(x),Q) is convex in F(x) so that Jensen’s inequality gives that 

k 

X~,s(F) < f inf /i(Q|P) + V'/if I f(x,s)ds\Apq(si)]dx. 

J A^QeVe ^ V 4 Si _! y 

Again by convexity of h and an application of Jensen’s inequality, we can further estimate 
f f(x,s)ds\Apq(s i ) S j < j h^f(x,s)\£pq(sj^ds. 

This proves X(F) < X(F). For the other direction we first consider the supremum over partitions 
S and let some S'-partition be fixed. The idea is to use a volume partition into equal sub-cubes 
with side length going to zero as a lower bound. More precisely, let {pfc(Z)}?= 1 denote the disjoint 
partition of A^ into cubes of volume \A' 1 \/k d and equal side length 5(k) = |A^_1 1 / d /k. Then, 

k d ^ r 

X~,s(F) > lim inf V -= inf h(Q |P) + V] h(F(p k (l) x (si-i,Si])\Apq(si)) 

k —^00 z ' K a QGVg z ' 1 

1=1 i=l 

= lim inf f inf h(Q\F) + h (fi (x)\Apq(si))dx 
k ^°° J Aj 

where each ff(x) is constant on each of the cubes pk(l), l = 1,..., k d . By Lebesgue’s theorem 
flf(x) —> J* 1 f(x,s)ds for Lebesgue almost all x as k tends to infinity. Therefore, by Fatou’s 
lemma and the lower semicontinuity of the rate function 

snpI~,s(F) > liminf [ inf h(Q\P) + h(fi(x)\Apq(si))dx 

H k=roc Q GVg f 

1 1=1 

> f liminf inf /i(Q|P)+ ^/i(/f(x)|A/r Q (si))dx 

J A' k =>-oo <Q> £Vg “ 

= / inf h(Q\F) + Y'/i( [ f(x, s)ds\Apq(si))dx = X S (F). 

J Aj QePfl ^ 
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For the supremum over S-partitions we use the same approach and consider a partition of 
intervals [0, k] for k G N with constant mesh size 1 jk. Using Fatou’s lemma, we have 


sup TgiF) > / lirninf inf 

S JA[ k^-oc QeV e 

ri/k 


k V 1 


i= 1 


i/k 


l/k 


dx, 


where we can look at k /( x, s)ds as a stepfunction on [0, k]. Similarly, for k\p^(i/k) - 


i - 1 )/k)} = k j^fJ.Q(s)ds with g ^ on [0, k], we have that 

-t pSj, rk 

^k h (. k / f{x,s)ds\kAg Q (i/k)^j = / ft(/*(r)|gj(r))dr. 


i=l 


'Si -1 


Now fix a: G Ai and fc > 1 and let Q k such that 

inf h(Q |P)+ /'"h(/ x fc (r)|4(r))dr = / l (Qf|P)+ ['h{f*(r)\gUr))dr, 

Q£7 J 0 JO JO 

which exists since lower-semicontinuous functions assume their minimum on compact sets. Let 
k n be the subsequence such that the limit inferior becomes a limit and for simplicity write again 
k. Further we can assume sup fc /i(Q^|P) < oo for Lebesgue almost all x since otherwise there 
is nothing to show. Since /i(-|P) has sequentially compact level sets there exists a cluster point 
Q* of (Qfc)fceN and by lower semicontinuity and Fatou’s lemma, we have 

p poo 

sup T S {F)> / h(Q* |P)+ / liminf/i(/^(r)|c/L(r))drda;. 

5 J a; Jo k 

Note that by Lebesgue’s theorem for almost all s G [0, oo), liminf^oo f£(s) = f(x,s). Further 
note that also liminf^oo g^ x («) = Indeed by the mean value theorem for s G 

((i — 1 )/k,i/k) there exists s' G ((* — 1 )/k,i/k) such that gqx.(s) = and 


s - 


dst 


>001 < 


d s^ 


:( s ') “ 


Win¬ 


ds' 


;00 “ ®/*QsOOI- 


The second summand on the right tends to zero as k tends to infinity by Corollary 17 For the 
first term we have by Lebesgue’s theorem 

~ sAmsj(s)I < ArQ^(E(^. _ i q(aI{o)~ 1 ) - g(s/(o) _1 )|X)) 


< Nw 2 |s' — s| 

which tends to zero as k tends to infinity and thus liminf^o Q gQx(s) = 
and lower semicontinuity gives 

lirninf h{f£(s)\gL{s)) > h(f(x, s)\£g Q ,{s)), 

/c —yoo K 

as required. 

Finally let F (J ACo- First, for any e > 0 there exists Qs,x G Vo such that 
1(F) > sup [ ^2 |Ah(s)|/i(Q5,x|P) 


;(s). Using this 


xdr. a 


F(A H (x)x(si_i,Si]) 


— e 


,s 

+ E|Ah(. 

1=1 

> sup [ i A =( ; 

x£S d i =1 

= s „ up [ S l As ( x )l sup[p F(Aa( [^ : ( ( 3 g|~ 1 ’ ai]) - (e p - l)A/iQ SiS (Si)] 

"i ,S . . . n 


F(A 3 (x)x(si_i,Si]) | a / 

|A= (x)| 


— £ 


— £ 


xeH d i=i 
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using also the Legendre transform of the relative entropy. Further, we have 

| VQs,*( s i) ~ ^ A R |Qs )X (r(sj,o) - r(s;_i,o)| < NX R W~ 2 \Si - S*_i|, 

and hence for p > 0 

r 

1(F) > p sup [ ^ ^[F(A h (x) x (si_i,sj]) - (e p - l)Ah\ R |AH(x)ll s i ~ «i-i|]] - £• 

x€3, d *=1 

If T is not right-continuous, there exists a point (x, s ) such that _F(x, s) < lim^^oo F(.T+l/n, s+ 
1/n) = M. Consider a sequence of finite partitions (E n , S' n )neN where the cube (n,=i ( x ji x h + 
1/n]) x (s, s + 1/n] is contained in (E n , 5 n ) for all n G N. Then 

1(F) > p[F(x + 1/n, s + 1/n) - F(x, s) - (e p - l)NX R l/n d+1 ] -e 

and letting n tend to infinity gives 1(F) > p[M — F(x, s)] — e which tends to infinity for p —> oo. 

If F is right-continuous but F ^ ACq there exists 5 > 0 and a sequence of measurable sets 
A k , with Vd+\(Ak) —> 0 and pp(A k ) > S. Using the regularity of the Lebesgue measure there 
exists a disjoint union of countably many d + 1-dimensional cuboids such that Ak C qf and 
v d +i((Ji Qi \ A k ) < Then, for every p > 0, 

OO 

1(F) > P J2 F (qi) ~ (e p - l)iVA R ^ +1 ((Jgf) - e 

i=i i 

> P^F(A k ) - (e p - l)NX R (n d+1 (A k ) + 1/k) - e. 

Letting k tend to infinity we have 1(F) > p5—e which tends to infinity as p tends to infinity. □ 


3.3. Contraction principle &: identification of rate function. In the present section, we 
apply the contraction principle to derive Theorem [2] from Proposition [9j Consider the function 
T : 1C — > M () given by 

F(-) ^ f({- x [0, oo)) fl {(y, s) € x [0,oo) : s < 

Then, the random measure \H(| AJ| _1 Y* ,t (- ) ■)) is exponentially equivalent to the random measure 
|A£| ~ l Y t (t~ 13 -). Moreover, T is continuous when restricted to the subset Li nc ,o(A / 1 x [0,oo)) of 
Lj nc (A / 1 x [0, oo)) consisting of those F with pp(dM) = 0, where 

M = {(y, s) <E A( x [0, oo) : s < \y\~ a } 

and 1C C Li nc (A / 1 x [0,oo)) denotes the family of all [0, oo)-valued, bounded, increasing and 
right-continuous functions. Since the Lebesgue measure of dM is 0, the rate function from 
the LDP of Proposition 9] is infinite on the complement of 1C. Hence, [TOl Lemma 4.1.5] shows 
that the random fields |ff Aj| _1 y* ,t (-, •)}t < i also satisfy an LDP on 1C. Hence, the contraction 
principle applies and it remains to identify the rate function. That is, we need to show that 

L^ h{m+ l K /fe - s)| s w(s, ) d5d! ' 

= f infh(Q|P) + h(g(y)\pQ(\y\~ a ))dy, 

Ja[ Q 

where / = dF/(dyds) and g = dG/dy denote the Radon-Nikodym derivatives of F and G, 
respectively. Note that if G was not absolutely continuous, then neither could be F, so that 
the left-hand side would be infinity. We show that the equality arises as a consequence of two 
inequalities. First, we consider the direction >. As in the proof of Proposition |9j an application 
of Jensen’s inequality implies that 

r\y\~ a u / r\y\~ a \ 

J h(f(y,s)\—HQ(s))ds > h(j f(y, s)ds\p Q (\y\~ a )j. 
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The right-hand side is equal to h(g(y)\pQ(\y\~ a )) if G(-) = 

The other direction is more involved. First, we proceed as in the proof of Proposition [9] and 
note that the right-hand side can be approximated using a suitable discretization. To be more 
precise, let {p(l)}f = 1 be a subdivision of A' x into congruent cubes of side length 5{k) = |A / 1 | 1 / rf 2 _fc . 
The point in the Z-th cube which minimizes the distance to the origin will be denoted by y k /. 
In the first step of the discretization, we replace the expression PQ(\y\~ a ) by p,Q(\y kt i\~ a ) ■ 


limsupY' [ inf h(Q\P) + h(g(y)\p,q(\y kt i\ a ))dy 

k —»oo Jp(l) 

< [ inf HQ\V) + h {9{y)\p®(\y\~ a ))dy. 

J A' QePe 

Proof. First, note that for every Z E {1, ..., 2 dk }, y E p(l) and Q E Tq we get that 

h{g(y)\m(\yk,i\~ a )) - h (9{y)\p®(\y\~ a )) 

< s(v) log +1 w(i»wr“) - Mtorn 

< \m(.\yk,i\~ a ) - m(\y\~ a )\ 

where the last inequality follows from the choice of y k .i. In particular, the right-hand side is 
always bounded above by 1. Moreover, for e > 0 we let A £ = {Z E {1,..., 2 dk } : min y^ p n) \y\ < 
e] denote the set of indices of cubes that are close to the origin. Then, the Lipschitz assumption 
implies that for every Z 0 A e , y E p(l) and <Q) E Vo, 

|mq(I? Jk,i\~ a ) - m(\y\~ a )\ < Na\y K i\~ a ~ l \y - y k f < (6) 

Hence, 

2 dk 

^(Ql 1 ?) + h(g{y)\p®{\y k ,i\~ a ))dy 

^ Jpd ) QePe 

~{J A , Q ig/(Q|P) + M^)|AiQ(|yr)) d y) 

< \p( 1 )\#A £ + Nae-^VdSik) < 2 d e d + JVa<r a_1 \/d£(A:), 

provided that k > 1 is sufficiently large. Since e > 0 was arbitrary, this completes the proof. □ 


The next lemma is proved similarly to Proposition [9] using Jensen’s inequality and a dis¬ 
cretization of the integral. We omit the proof. 

Lemma 19. Let k > 1 and 1 < l < 2 dk be arbitrary. Then, 

IpWT 1 [ „ in | MQI p ) + h{g(y)\p®(\y k ,i\~ a ))dy 

Jp(l) 

> inf h(Q |P) + / l (|p(Z)|- 1 G(p(Z))|^(|y fc) r“))dy. 


Now that we have discretized the integral, we can define approximations to the desired 
function F. For this purpose, we first need to construct certain minimizers. Recall from 

l ) is continuous on every set of the form {Q : 


Corollary 17 that the function Q i-)- fiQ(\y k ,i\ 

< K\ for some K < oo. Therefore, the function 


Q h(Q |P) + h(|p(Z)r 1 G(p(Z))| / r Q (|y fc , / |-“)) 
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is lower semicontinuous, and we let Q k ,l be one of its minimizers. Now, define measurable 
functions /^ : A\ -> [0, oo], k > 1 by 


f {k) (y,s) 


|p(f)|~ 1 G(p(I)) ^AXQfc.i («) 

VQ ktl (,\yk,i\- a ) 

.T s m k A s ) 


if V e p(l) and s < \y k ,i\ “, 
if y G p(l) and s > | y k ,i\~ a - 


Here, we make the convention that the first line is equal to zero if pq k , (\y k ,i\~ a ) = G(p(l )) = 0 
and is equal to infinity if pq k ;d2/fc,zl _a ) = 0, but G(p(l)) 0. Furthermore, we let F^ denotes 

the distribution function of the measure with density f^ k \y , s ). Then, for every y £ p(l), 


inf 

Q6P 9 


/i(Q|P)+ f h^f {k \y,s)\^pQ(s)^jds 


< h(Q k ,i\¥) + 






1 d-s 


,(s))ds 


= h(Q k>l \F) + h(\p(l)\- l G{p{l))\pq Kl (\y K z |"“)). 

By the goodness of the rate function in Proposition [ 9 J the functions (F^)k> 1 have an accu¬ 
mulation point, and lower-semicontinuity therefore gives that 


/ inf 
/a; 


+ / h(f{y,s)\p®(s))ds < / inf h(Q\F) + h(g(y)\pq(\y\ “))d y. 

Jo Ja[ Qe:P e 


Hence, it remains to show that the measure induced by G(-) coincides with the measure induced 
by F(-1m)- In order to prove this claim, we first show that F ( ' k \M lyk ^ \ M ) tends to zero as k 
tends to infinity, where M^> = {(y,s) € A\ x [0, 00 ) : y € p(l ) and s < \y k ,iA a } 


Lemma 20. The expression F^ k \M^ \ M ) tends to zero as k tends to infinity. 


Proof. First, observe that F^ k \M^ \ M ) can be expressed as 


F^{M^\M) = Y J G{p(l))W)\ 


-1 




1=1 


'p(i) 


MQfc-d yk,l\ “) 


-d y. 


Now, for e > 0 introduce the set 

A e = {l £ {1,..., 2 dk } : min \y\ < e or max \y\~ a > tcs m i n - e} 

y&p{l) y&p(i) 


of indices whose associated cubes are far away from the origin and the boundary of the ball 
B,„ a \-i/ a (o). Hence, we arrive at 

F (fc)( M (fc) \ M ) = ae~ a ~ l NVd5{k) V -+ r ( e ) ( 7 ) 

MQ M d2/fc,zl “) 

where r(e) tends to zero as e tends to zero. Since the sum above consists of at most 2 dk 
summands, it suffice to consider those l 0 A e that satisfy \p(l)\~ 1 G(p(l)) > pq k t (\y k ,i\~ a )- 
Now, note that if l 0 A s , then p^{\y k ,i\~ a ) > 1/IF for some sufficiently large K = K{e) not 
depending on k,l. We claim that also pq k l (\y k ,i\~ a ) > 1/IF. Once this is shown, the proof is 
complete. 

Suppose that pq k ^(\yk,l\~ a ) < l / K - First > if \pQ')\~ 1 G(p(l)) > PF(\yk,l\~ a ), then 

h(\p(l)\~ lG (p(l))\p ¥ (\y k fi- a )) < h(\p(l)\- 1 G(p(l))\pq kl (\y k ,i\~ a ))^ 

which contradicts the minimality of Q k ,l- Otherwise put Q* = AP + (1 — A )Q k ,l> where A £ [0,1) 
is chosen such that 


^M\yk,i\ a ) + (i-A) MQfe XM a ) = \p(i)\ 1 g(p(i)). 
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Then, since the specific relative entropy h (introduced in Section [TJ is an affine function, 

kq * ip) + / l (ip(or 1 G( /3 (/))| MQ .(| 2/fc , z |-“)) 

= (l-A)h(Q fc)/ |P) 

< h(® k>l |P) + h{\p(i)\- 1 G( P mm k A\yk,ir a ))> 

which contradicts again the minimality of Qk,i- □ 


Now, we can complete the proof of Theorem [2] 

Proof of Theorem ^ By Lemmas 18 and 19, it suffices to show that G(f) = F(f 1m) holds for 
any / : A] —> [0, oo) of the form / = l p (; 0 ) for some 1 < Iq < 2 dk °. Now, we can argue as follows. 

\F(fl M ) - G(f) | < limsup|F«(/l M(fe) ) - G(/)| +hmsu Pj F^(/l M(fe)w ). 


k —^oo 


k—>oo 


Lemma 20 shows that the second summand is zero. Moreover, by definition of F^ k \ 

F ( - k Hfl M w) = GW(f), 

where G^ k \ k > fco denotes the measure with locally constant density g^ k \y) = \p(l)\~ 1 G(p(l)). 
Hence, 

GW(f) = GW(p(lo))= J2 G{p(l)) = G(p(l 0 )) = G(f), 

p(l)cp(lo) 

as required. □ 


3.4. Proof of Corollary [3[ 


Proof. The upper estimate is a direct consequence of the upper bound in Theorem [2] 

lirnsup |A(| _1 logpt = limsup |A£| _1 logP(|Aj| _1 y t (i _|S -) = 0) 
t-> o t->o 

< - / mf (h(Q |P) + X R Q(T(\y\- a , o)))dy. 

J A' QePe 


For the lower estimate first note that 


Pt = E exp 



where T(a,y) = E (q(al(y) 1 )|^f) is a non-local function of the transmitter process. In order to 
be able to apply m Theorem 3.1], we need to establish a translation-invariant setting using 
discretization of the integrand. To be more precise, we subdivide A' t into 2 dn sub-cubes A\ of 
side length 2 1 ^ n |u;s m i n t| _/3 and let yt denote the corresponding element of the subcube A \ which 
is closest to the origin. Then 


Pt> E exp ( — A 


2 dn 

£ 

i=1 • 


2 dn . 

r(i _1 ^(|yi|),y)d7/) >Eexp(-A R ]E / T^r 1 C(\yi\) 

i =l •' K t 


where T b (a,y) = E (q(al b (y) 1 )|A"). Further let A" per ’* be the conhguration obtained after 
extending the configuration of the marked Poisson point process X in the subcube A\ peri¬ 
odically in the entire Euclidean space M. d . The error made replacing X^=i Ja* y) 

by Ei=i /a ! T per ’ ?l (t _1 f(2/i), y) is negligible in the large deviation principle where r per,fe (a,y) = 
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E(g(a/ b (y)- 1 )|XP er - i ), indeed 

f [T b (t-H( yi ),y) - rP er ’ fe (i- 1 f(y i ),y)]dy 
J hi 


< N(w*t)- L / [/>, y, X ) - />, y, X^))dy 
J hi 


< N(s 2 mi „w 2 t) 1 


E 

x,ex 

Xj^h l t ,Xjdh l t ^ b 


£ b (\Xj - y\)dy + / 4(|Xj - y|)dy 


x j e.YP er ’ ! 

xrfhlxjehi. 


< AT(4 in te4)-MX(A* i6 \ AJ) + X^\A\ >b \ A*)] / 4(M)dy, 


B b (o) 


( 8 ) 


where AJ b denotes the volume A\ joined with its 6-boundary. Consequently, for all e > 0, 

2 dn 

Eexp(-A R ^ f T b (t~ l l{\ yi \ 

i= 1 

2 dn 

> e-^'Eexp (- a r ^ f r^rto^i),: 

»=i 

2 dn 

- p (a R E /, r b r ^( 2 /i), y) - r^(rt( W ),y)di, > e| a;i) . 

j=l ■' A t 

By equation Q, the second line is bounded from below by 

-2P(A R 21V( S 2 min a; 2 t)- 1 f £ b (\y\)dyX (Aj b \ Aj) > e|Aj|). 

V JB b {o) 7 

But this goes to zero on an exponential scale infinitely fast by Lemma [4] Hence, 


dm |Aj| 1 logEexp ( - ArV / T b (t 1 £(y i ), 

*-° V ->K 

2 dn „ 

> lim|A'|- 1 logEexp(-A R y; / T^ b (r 1 £( yi ), 

2 dn 

= 2- d ^ lim |Ar 1 logEexp(-A R [ 

z —< t-* o v \i 


where we used the independence of the X per ’’ 1 with respect to i in the second line. Now we are 
in the position to apply m Theorem 3 . 1 ] and write 

2 dn 

jim |Aj| _1 logEexp ( - Ar ^ / r per ’ b (rt(y;), y)dy ) 

’ i= 1 ^ A i 

2 dn 

> —2~ dn mf (A(Q|P) + A R Q(r 6 (|yj| _Q , o))), 

2=1 

where we also used the continuity of r per,b ensured by the truncation of the interference. Notice 
that 


limsup inf (h(Q|P) + A R Q(r b (|yj| _ ", o))) < inf (h(Q|P) + A R Q(r(|yj| -Q , o))). 

6_>oo QeVo 
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Indeed, let Qo be a minimizer of the right hand side, then 

limsup inf (/i(Q|P) + \ R Q(T b (\yi\~ a , o))) < /i(Q 0 |P) + A R limsupQ 0 (r i (|y i | _a , o)) 

oo 

and it suffices to show that 

|Q 0 (r 6 (|y|-“,o))-Qo(r(|2/|-“, 0 ))| <IV|y|-“ U ;- 2 Q 0 ( £ ^(|^|)) 

Xjt£A b 

tends to zero as b tends to infinity. But this is true since Qo is a translation-invariant point 
process. In order to perform the large-n limit, we have to show that 

V* inf (h(Q\¥) + X R Q(T(\y\- a ,o))) 
is continuous. But this is also true since 


inf {h(Q\F) + X R Q(T(\y\- a ,o)))- inf (h(Q\F) + A R Q(r(|x|-“, o))) \ 

QeVg QePe 

< sup | A R Q(r(|y| _ ", o)) — A R Q(r(|a:| _a , o))| < iVu; _1 A R ||y| _a — \x\~ a 

Qe-Po 


This gives the result. 


□ 


4. Importance Sampling 

In this section, we show how the LDPs derived in Theorems [l] and [2] can be used to devise 
an importance-sampling scheme improving the accuracy of basic Monte Carlo approaches for 
estimating the probability of observing unlikely configurations of connectable receivers. Theo¬ 
rems [l] and [2] imply that such probabilities generally tend to zero exponentially quickly, so that 
basic Monte Carlo estimators perform poorly. 

The general heuristic for devising importance-sampling schemes is the following. Instead of 
sampling the transmitters according to their true distribution, the simulation is performed by 
using a modified law under which the considered rare event is more likely. An appropriate 
reweighting using likelihood ratios ensures the unbiasedness of the new estimator. For a more 
detailed discussion of the general technique of importance sampling, we refer to the textbooks f2j 

m- 

In principle, Theorems [l] and [2] provide precise descriptions of the asymptotically exponen¬ 
tially optimal change of measure, in the sense that the modified law of transmitters should be 
given by suitable Gibbs point processes. However, as these distributions just arise as minimizers 
of fairly complicated functionals, it is difficult to use them for computational purposes. Still, by 
performing this minimization in the restricted class of Poisson point processes, we can achieve 
substantial accuracy benefits. 

We only provide a proof-of-concept for the use of importance sampling, and therefore assume 
a specific parameter constellation in the following. First, we fix d = 2, w = A R = Xj = 1 
and assume that the path-loss function is given by l(r) = r -4 . Moreover, we assume that 
there is no random environment, and that transmission powers and fading random variables 
are constant and equal to 1. Note that this choice is not covered by the assumptions for The¬ 
orems [l] and [2j Nevertheless, our simulation results illustrate that variance reduction through 
importance sampling also hold under weaker conditions than the ones assumed in Theorems [T] 
and[2j 

4.1. Importance sampling related to Theorem [lj Since we have assumed that there is 
no random environment and that transmission powers and fading variables are constant, the 
minimization in the rate function of Theorem [l] is performed only over stationary point processes 
of transmitters and receivers. As mentioned above, this minimization is intractable in its full 
generality. Nevertheless, in this section, we show that if minimization is performed only in 
the class of Poisson point processes, then the problem becomes tractable. In fact, we provide 
an example problem, where the minimization can be reduced to a standard two-dimensional 
constrained minimization problem, where the constraint is given in terms of certain special 
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functions. The disadvantage of this approach is that solving the minimization problem in a 
restricted class of point process will not automatically lead to good choices for the importance 
sampling. This will become apparent from the simulation results discussed below. 

We assume that t = 1 and consider events of the form 

An ' a = {\h\ ^2 < a }, 

Xi£ A n 

i.e., the event that the (spatially) averaged number of connectable receivers associated with 
transmitters in the cube A n is less than a. Now, we explain how to implement an importance¬ 
sampling scheme based on the LDP. A related importance-sampling scheme for a Poisson point 
process on the real line has already been considered in but for the convenience of the 
reader, we present some details in our situation. 

In order to estimate the probability of the event A n a , we simulate the Poisson point processes 
with new intensities pr > 0 and pj > 0 in A n . Then, the likelihood ratio of Poisson point 
processes with intensity 1 with respect to these point processes is given by 

exp (|A n |(/i R - 1) + |A n |(/x T - 1 ))p R X(An 'W 1 (An) - 
Hence, an unbiased estimator for P(A nja ) is given by 

Pn,a,n r ,» R = exp (lAnK/ZR - 1) + |A n |(jU T ~ 1)) Mr^ ^^ 

Note that in order to take into account edge effects, we also generate transmitters and receivers 
with the unmodified intensity in a small environment around A n . We can obtain estimates p and 
v of the expectation and variance of p n ,a,fi T ,fj . R by considering the sample average and variance 
of IV > 1 independent copies generated using Monte Carlo simulation. 

This leaves the question as of how to find good choices for pr and pj. In a first attempt, 
choose these parameters to minimize the large-deviation rate function appearing in Theorem [lj 
If Q E Vg the distribution of independent Poisson point processes of receivers and transmitters 
with intensities pr and pj, then the relative entropy /i(Q|P) is given by the formula 


{PR log PR - MR + 1) + (hT log /XT - PT + I)- (9) 

Hence, to determine the optimal intensities (AR. opt , Ay, opt) according to Theorem [TJ we need to 
minimize Q under the constraint 

Q*(#y(°)) < a. (10) 


Next, we express the constraint (10) in terms of certain special functions. First, by Campbell’s 
theorem, 

I —4 


"(#^ ( °)) = PtPr 


> 1) da; 


Bi(o) + X)i>l l-^il 4 

PtPr2tt f rp( V] < p\(r A - l))dr, 

• /o V i>i 7 


where in the last line we used that scaling by 1 /y/pr transforms a Poisson point process with 
intensity 1 to a Poisson point process with intensity pj. Moreover, X)i>i l^*l~ 4 i s distributed 
according to an inverse gamma distribution with parameters 0.5 and 7r 3 /4. In particular, 
P(Ei>i < «) = ^“ 1/2 t(i/2 , —7t 3 /(4s)), where denotes the incomplete gamma 

function. Now it is easy to check that 

Q*(#Y (o) ) = PRPT27T f rj(l/2, 7r 3 / (Apj(r 4 - l)))dr 

Jo 

= prpjtt exp (7r 3 /(4^T))erfc(vr 3/2 /(2^T)), 

where erfc denotes the complimentary error function. For instance, if we choose a = 0.5, then 
(Ar, opt, Aj, opt) ~ (0.832,0.984). 
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(/^Rj P-t ) P V 

(1,1) 3.31 x 10 -4 3.31 x 10 -4 

(0.832,0.984) 3.12 x 10" 4 3.69 x 10" 4 

(0.892,0.989) 3.29 x 10" 4 5.10 x 10" 5 

Table 1. Comparison of the simulation results for the expectation and variance 
of the considered importance sampling estimators with transmitter and receiver 
intensities (p,R,pj). 


In order to assess the actual accuracy improvements that can be achieved with this importance 
sampling-scheme, we performed a prototypical Monte Carlo analysis. We fixed a = 0.5, n = 25 
and performed N = 1,000,000 simulation runs. We consider three different parameter choices 
for the importance sampling intensities (/zr, /zj). First, we consider the case of basic Monte Carlo 
simulation, that is (/z r, pj) = (1,1). Second, we take the intensities that are obtained from the 
large-deviation analysis performed above, i.e., (/z r, pj) = (0.832,0.984). Third, we estimate 
(/zr,/zj) from a simple cross-entropy scheme. That is, we performed a pilot run of 100,000 
basic Monte Carlo simulations and determined the average intensities under the condition that 
the rare event occurs. This gives (/zr,/zj) = (0.892,0.989). We refer the reader to (19j for 
further details on the general cross-entropy technique. The results for p and v are reported in 
Tabled! 

In particular, we would like to draw the attention to an important observation: The estimator 
that is obtained as the solution of the optimization based on our large-deviation principle 
actually has a higher variance than the basic Monte Carlo estimator. Given the close relation 
between large-deviation theory and asymptotically optimal change of measures, this might come 
as a surprise at first sight. However, since we performed our optimization not in the full class of 
stationary point processes, but only considered Poisson point processes, the simulation output 
does not contradict this intuition. In fact, considered from a different perspective, the simulation 
results provide evidence that the optimal change of measure is rather far (in the Kullback-Leibler 
distance) from being a Poisson point process. In contrast, performing the change of measure 
with the intensities obtained from the pilot run shows that for the considered example, a more 
than seven-fold variance reduction can be achieved. 

The discussion in the previous paragraph raises the legitimate question whether the change 
of measures deduced from the large-deviation result are of any practical use for importance 
sampling? Indeed, in the example described above, the intensities that lead to the seven-fold 
decrease in variance could be found without reference to the LDP, namely by an ‘educated 
guess’ (or rather ‘cross-entropy’). Nevertheless, when considering importance sampling in the 
setting of Corollary [3] finding a good importance sampling change of measure would involve 
‘guessing’ a continuous family of parameters, which is substantially more involved than what 
we have done above. In contrast, a simple analysis of the large-deviation rate function provides 
immediately a useful heuristic for the shape of the curve. 


4.2. Importance sampling related to Theorem [2], Finally, we investigate importance- 
sampling techniques related to Theorem [2j We consider the specific setting of Corollary [3j i.e., 
estimation of the isolation probability pt for small values of t. Similar to the situation considered 
in Section 4.1, the full minimization problem is intractable, so that we restrict our attention to 
the class of homogeneous Poisson point processes. However, the situation is slightly different 
from the one considered in Section 4.1 Instead of globally optimizing a transmitter and receiver 
intensity, we now have the freedom to choose a different intensity for each point in A). Due 
to isotropy, this reduces to the task of choosing an optimal intensity A opt (r) for each r e [0,1]. 
This optimal intensity must minimize the following expression that can be derived from the 
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V 

A(-) = 1 7.72 x 10 _b 8.22 x 10~ iU 

A(-) = Aopt(-) 7.70 x 10~ 6 1.78 x 10" 10 
Table 2. Comparison of the simulation results for the expectation and variance 
of the basic versus the importance sampling estimator. 


variational characterization in Corollary [3| 

A 0 pt(r)logA opt (r) - A opt (r) + 1 + P(r -4 > 1 + A opt (r) 2 ^ |W|" 4 ). ( n ) 

i> 1 

This is a standard minimization problem that can be solved by finding the roots of the derivative 
with respect to A op t(r). After some simplifications, we arrive at 

, x XX tt ,7T 3 A op t (r) 2 

logAopt(r) = > /F^I eXP( 4(r- 4 -l) ) - 

This equation can be solved numerically; a plot of this solution is shown in Figure [2j 



Figure 2. Plot of the optimal density A opt (r) at distance r from the origin 

As in the previous example, in order to assess the actual accuracy improvements for the 
estimation of pt that can be achieved with this importance sampling-scheme, we performed a 
prototypical Monte Carlo analysis. We fixed t = 0.002 and performed N = 1,000, 000 simulation 
runs. 

In contrast to the previous example, we see that the importance-sampling estimator derived 
from large-deviation theory provides substantial benefits. Indeed, the variance is reduced by 
approximately 78%. Furthermore, applying the cross-entropy technique for the present example 
would be substantially more involved than in the previous one. Indeed, instead of simply 
estimating two parameters, we would need to extract an entire curve from the pilot runs, so 
that proper statistical tools would be needed to estimate such a functional object from data. 
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